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ABSTRACT 



12 We have developed a general model for determining density-dependent effec- 

tive dielectronic recombination (DR) rate coefficients in order to explore finite- 
density effects on the ionization balance of plasmas. Our model consists of mul- 
tiplying by a suppression factor those highly-accurate total zero-density DR rate 
coefficients which have been produced from state-of-the-art theoretical calcula- 
tions and which have been benchmarked by experiment. The suppression factor 
is based-upon earlier detailed collision-radiative calculations which were made for 
a wide range of ions at various densities and temperatures, but used a simplified 
treatment of DR. A general suppression formula is then developed as a func- 
tion of isoelectronic sequence, charge, density, and temperature. These density- 
dependent effective DR rate coefficients are then used in the plasma simulation 
code Cloudy to compute ionization balance curves for both coUisionally ionized 
and photoionized plasmas at very low (rie = 1 cm~^) and finite {uc = 10^^ cm~^) 
densities. We find that the denser case is significantly more ionized due to sup- 
pression of DR, warranting further studies of density effects on DR by detailed 
coUisional-radiative calculations which utilize state-of-the-art partial DR rate co- 
efficients. This is expected to impact the predictions of the ionization balance in 
denser cosmic gases such as those found in nova and supernova shells, accretion 
disks, and the broad emission line regions in active galactic nuclei. 



13 Subject headings: suppression 
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Introduction 



15 Astronomical emission or absorption sources have an enormous range of densities. 

16 Two examples include the intergalactic medium, with ~ 10~^ cm~^, and the broad 

17 emission-line regions of Active Galactic Nuclei, with Up ^ 10^° cm~^. The gas pr oducing 



IB the spectrum is not in thermodynamic equilibrium ( lOsterbrock fc Ferland 



20061 ). so 



19 microphysical processes determine the physical conditions. 

20 The two common cases encountered for ionization are photoionization and collisional 

21 (e.g., electron-impact) ionization. In both cases, ions are recombined by dielectronic 

22 and radiative recombination, with dielectronic recombination (DR) usually the dominant 

23 process for elements heavier than helium. Databases give ionization and r ecombination 

24 rates that are the sum of sever al contributing proc e sses. Examples include 



25 for electron impact ionization, 



26 DR project (IBadnell et al. 



Voronovl (119971 ) 



Verner fc YakovlevI (119951 ) for photoi onization, and t he 



2003) for dielectronic recombination and Badnell ( 2006a ) for 





27 radiative recombination; it is these latter data □ which will be of primary interest to us in 

28 the present study. 

29 The collisional ionization and recombination rate coefficients used in astrophysics are 

30 frequently assumed to depend on temperature but to have no density dependence. The 

31 rigorous treatment of density dependent ionization and recombination rate c oefficients is 



32 via collisional-radiative modeling. This was introduced by 



Bates et al. 



( 119621 ) for radiative 



33 recombination onl y and extended to tr e at th e much more complex case of dielectronic 



34 recombination by 



Burgess fc Summers! ( 119691 ). Summers applied their techniques to 



35 determine density dependent ionization and recombination rate coefficients, and the 

36 consequential ionization balance for collisional plasmas, for H-like thru Ar-like ions. 



^http : / /aindpp . phys . strath . ac . uk/tamoc/DATA/ 
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37 Graphical resu lts were present ed for the elements C, O and Ne ( 1 Summers! 1 1 9 72l ) and then 



38 N, Mg and Si ([Summers 



19741 ) ■ Reduced temperatures and densities were used so as to 
enable easy interpolation for other elements. Tables o f such reco r nbination rat e coefficients 
were made available only via a Laboratory Report — ISummerd (11974 fc 19791 ) — due to 



their voluminous nature at that point in history. The 'difficulty' in utilizing this pioneering 
data le d to so me modelers attempting to develop simplified approaches. For example, 



Jordan! (11969! ) used an approach which was based on truncating the zero-density DR sum 



over Rydberg states using a simple dens ity dependent cut-off which itself was based on 



45 early collisional-radiative calculations by 



Burgess fc Summers! (11969! ): a suppression factor 



46 was formed from its ratio to the zero-density value and then u sed more generally. Also, 



Davidson! ( !l975l ) simplified the collisional-radiative approach of iBurgess fc Summers! (11969! ) 



48 and, using hydrogenic atomic data, determin ed suppression 



"actors for Li-like C IV and 



49 O VL New calculations for C IV were made by; 



(generalized) collisional-radiative modeling (jSummers Sz Hooper 



Badnell et al 



(!l993!l ut ilizing more advanced 



atomic data at coUisional plasma temperatures (see the references in 



1983[) and much improve d 



Badnell et al. 



(119931)). 



52 All of the above works were 

53 DR data (excluding 



Badnell et al. 



o r elec tron collisional plasma s and us e d rat 



ler basic 



( 11993! )) as epitomized in the iBurgessI (!l965! ) General 

54 Formula, viz. a common dipole transition for dielectronic capture, autoionization, and 

55 radiative stabilization. The purpose of the present paper is to explore density suppression 

56 of DR in photoionized plasmas, and within collisional plasmas, using state-of-the art DR 

57 data which takes account of a myriad of pathways not feasible in the early works above, 

58 but which has been shown to be necessary by comparison with experiment. We wish to 

59 gain a broad overview utilizing the large test-suite maintained by the plasma simulation 



60 code Cloudy. We utilize an approach to DR 

61 detailed collisional-radiative results given in 



suppression w h ich is motivated initially by the 



Badnell et al 



(119931) for C IV at T = 10^ K, 



62 along with known scalings to all temperatures, charges, and densities. Using these results as 
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63 a guideline, a more general suppression formula is then determin ed by fitting to suppres sion 



1974 fc 1979h for 



64 results from extensive detailed coUisional-radiative calculations ( ISummersI 

65 a wide range of ions at several densities and (high) temperatures. Additional modifications 

66 are then introduced to account for low temperature DR. 

67 The outline of the rest of the paper is as follows: in the next section we describe the 
6B DR suppression model we use; we then apply this suppression to the zero-density DR data, 

69 and use the resultant density-dependent DR data in Cloudy to determine the ionization 

70 distribution produced under photoionized and collisional ionization equilibrium at low and 

71 moderate densities. 



2. Generalized Density Suppression Model 



We use the following approach, detailed more fully in th e subsections below. First 



Badnell et al 



f ll993h for DR 



74 the high-temperature coUisional-radiative modeling results of 

75 suppression in C IV are parameterized by a pseudo-Voigt profile to study the qualitative 

76 behavior of suppression as a function of density and temperature. Next, this formulation 

77 is then used as a guideline for developing a more comprehensive suppression formula which 
7B is obtained by fitting to collisional ra diative data for various isoelectronic sequences, ionic 



79 charges, densities, and temperatures ( Summer; 



1974 fc 19791 ). Lastly, the suppression 



80 formulation is extended to low-temperatures according to the nature of the sequence-specific 

81 DR. 



82 2.1. High- Temperature Suppression for Li-like C IV 

83 We begin by considering DR of Li-like C IV, for which the density dependent total DR 

84 rate coefficient, and therefore the suppression factor, has been computed rigorously within 
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85 a collisional-radiative modeling approach f lBadnell et al.lll993l ) 



86 In the electron coUisional ionization case, because of the consequential high temperature 

87 of peak abundance, dielectronic recombination occurs mainly through energetically high- 

88 lying autoionizing states (via dipole core-excitations) for which radiative stabilization is by 

89 the core electron into final states just below the ionization limit: 

e" + ls^2s ls^2pnl -> ls^2snl + hu . (1) 

90 In the zero-density limit, the intermediate ls^2sn/ states can only decay further via radiative 

91 cascading until the ls^2s^ final recombined ground state is reached, thereby completing the 

92 DR process: 



ls^2snl 



ls^2sn'l' + hvi ... ls^2s^ + hvi + hi'2 + ... 



(2) 



93 For finite electron densities rze, on the other hand, there is also the possibility for reionization 

94 via electron impact, either directly or stepwise. 



e + ls^2sn/ 



ls^2sn'l' + e- ^ ... ^ ls^2s + e" + e" , 



(3) 



95 and the probability of the latter pathway is proportional to the electron density n^. Because 

96 of this alternative reionization pathway at finite densities, the effective DR rate coefficient 

97 a5^(ne, T) is thus suppressed from the zero-density value aD^^{T) by a density- dependent 

98 suppression factor 5 (tt-c, T): 



aDR(ne,T) = S'(nc,T)aDR(T) . 



(4) 



From the earlier detailed studies of 



Davidson! (119751 ) and 



Badnell et al. 



(119931 ). the 



100 suppression factor is found to remain unity, corresponding to zero suppression, at lower 

101 densities until a certain activation density rig a is reached, beyond which this factor decreases 

102 exponentially from unity with increasing density. We have found that this suppression 
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factor, as a function of the dimensionless lo g density parameter x = log^^Q rzc, can be modeled 



104 quite effectively by a pseudo-Voigt profile (IWertheim et al. 



19741 ) — a weighted mixture /i 



105 of Lorentzian and Gaussian profiles of widths w for densities above the activation density 



Xa = login nc.a: 



S (x; Xa 



1 



l+(£^)2 

Fitting this expression to the suppression factor of 



X <Xa. 
X > X^, 



(5) 



Badnell et ali fll993h for C IV (which 
was computed at T = 10^ K) yielded the values /i = 0.372, w = 4.969, and x^ = 0.608, 
and this parameterization formula is found to be accurate to within 5% for all densities 
considered (see Fig. [1]). 



2.2. Generalized High- Temperature Suppression Formula 



Given the suppression formula for Li-like C IV, corresponding to ionic charge Qq = 3 
and temperature Tq = 10^ K, we wish to generalize this expression to other Li-like ions of 



charge q and (high) T according to the follow i ng qu alitat ive guidelines. It is we 



density effects scale as q'^ — see 



Bates et al. 



1 known that 



(119621 ) and lBurgess &: SummersI (119691 ). The 
activation density is attained when the reionization rate in Eq. |3l which depends linearly on 
the density, becomes comparable to the radiative stabilization rate in Eq. O The radiative 
rate is independent of density and temperature, but scales with charge as A,. ~ g^, whereas 
the electron-impact ionization rate depends on all three, viz. rieaeii ~ ncq~^T~^/'^ . An 
initial suggestion is that the activation density is attained when these two are approximately 
equal, i.e., 

ne,ag-^T-V2 ^ , (6) 
indicating that the activation density should scale as riea ~ g'^T^/^, if the above qualitative 



123 discussion holds. The log activation density for all q and T might therefore be expected to 
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Fig. 1. — Pseudo-Voigt fit of tfie suppression factor for C IV, as given in Eq. |5]with a scaled 
activation density as given by Eq. [TJ shown for two different temperatures. The red solid 
curve shows that the parameterization for T = 1 x 10^ K, corresponding to an activation 
density of x^. = 0.608 ( with = 0.372 and w = 4.969), is in close agreement with the actual 



data of 



Badnell et al 



(119931 ) (open circles). The blue dashed curve is the parameterization 
for T = 1.5 X 10^ K, using instead an activation densit y of x„ = 0.196 ( and the same /x and 



w), and giving satisfactory agreement with the data of 



Davidson! (jl975[ ) (solid circles) 
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124 obey the scaling relationship 



(go, To) + logio 



T 



1/2 



(7) 



125 where Xa{qo, Tq) = 0.608, go = 3, and To = 10^ K a re the (log) act i vatio n density, the charge 



126 and the temperature for the C IV case treated by 



Badnell et al. 



(119931 ). We note that this 



127 expression, when applied to Li-like O VI, gives an increase in the activation den sity by a 



12B factor of (5/3)'' = 35.7, in agreement with the approximate factor of 40 found by 



Davidson 



129 f ll975l ). Furthermore, when sea 

130 with the suppression results of 



ed in temperatu re, the formula gives fairly good agreement 



DavidsonI f ll975h for C IV at T = 1.5 x 10^ K (see Fig. [I]). 



2.2.1. Fit to the Collisional Radiative Data 



The preceding treatment reasonably extends the C IV suppression factor at 10^ K 
to other high temperatures and to other Li-like ions. However, we need suppression 
factors applicable to all ionization stages of all elements up to at least Zn for a general 
implementation within Cloudy. Unfortunately, detailed coUisional-radiative modeling data 
with state-of-the-art DR data is still rather limited. H owever, ex t ensive tables of effective 



SummersI (11974 &: 19791 ) for a wide 



recombination rate coefficients have been computed by 
variety of isoelectronic sequences, charge-states, temperatures, and densities. The treatment 
of DR there is somewhat simplified, but we only require the ratio of finite- to zero-density 
rate coefficients to determine the suppression factor. We then combine this ratio with our 
state-of-the-art zero density DR rate coefficients again for use within Cloudy. This ratio is 
much less sensitive to the specific treatment of DR. 

The rather simplistic scaling formula in Eq. [7] was found to be inadequate when 



144 applied to the extensive tabulation of suppression factors found in 



SummersI ( 1l974 fc 19791 ). 



145 Instead, in order to fit the data accurately, a more generalized formula was arrived at. 
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146 where a pseudo-Gaussian, corresponding to /i = in Eq. [5], was more appropriate, 



S^{x;q,T) 



1 



X < Xa{q,T, N) 



, x-Xa{q,T,N) n2 



(8) 



Furthermore, the activation density was found to be best represented by the function 

7 / rn \ 1/2" 



Xa{q,T,N) 



X" 



log 



10 



T 



qo{q,N)J yTo{q,N) 

where the variables qQ{q,N) and To(g, A^) are taken to be functions of the ch arge q and 
the isoelectro nic sequence, labeled by A^. A fit of the suppression factors of 



(9) 



Summers 



(11974 fc 19791 ) for all ions yielded a global (log) activation density x^ = 10.1821 and more 
complicated expressions for the zero-point temperature Tq and charge go- These were found 
to depend on both the ionic charge q and the isoelectronic sequence A^ viz. 



To{q,N) = 5 X 10^ [qo{q,N)f 



(10) 



and 



qo{q,N) 



;i - VWq)A{N)/^ 



'IV. 



where 



A{N) = 12 + lOA^i + ^^^^^i — '^(N-Ni) 



(12) 



N1-N2 

155 depends on the isoelectronic sequence in the periodic table according to the specification of 

156 the parameters 

/(3, 10) A^G2"'^row (37,54) A^GS^'^row^ 
(A^i,A^2) = (11,18) A^GS^'^row (55,86) A^GG'^^row • (13) 
\(19, 36) A^ G 4*^ row (87, 118) A^ G 7*'' row j 

157 However, even this rather complicated parameterization was inadequate for the lower 

158 isoelectronic sequences A^ < 5, and for these we explicitly list the optimal values for A{N) 
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159 in Table [TJ Furthermore, at electron temperatures and/or ionic charges for which the 

160 g-scaled temperature 9 = T /q^ was very low {9 < 2.5 x 10^ K), a further modification to the 

161 coefficients A{N) for < 5 is necessary in that the values in Table [T] should be multiplied 

162 by a factor of two. 

163 The above final formulation, which consists of the use of Eq. [HI with /i = 0, w = 5.64548, 

164 and a rather complicated activation density given by Eqs. |9l [TOl [TTl [121 and [131 with 



165 X 







10.1821, has been fo und to mo d el the entire database of ions, temperatures, and 



166 densities considered in the 



Summers! (11974 fc 19791 ) data fairly well. To illustrate the 



167 general level of agreement over a large range of ions and environments, we comp are our 



168 parameterize d model formulation to the actual suppression data from that report ([Summers 



1974 fc 1979[ ) for a few selected cases in Fig. [2l In order to quantify more fully the extent of 



agreement, we focus on the case of iron ions, for which we study density effects on ionization 
balance determination in the next section. A comp arison is shown in Fig. [3l between our 



172 predicted suppression factors and the data from the 



Summers! (11974 fc 1979[ ) report. It is 



173 seen that our m odel fits that data to wi thin 21% for all densities, temperatures, and ionic 



stages reported ([Summers 



1974 fc 1979[ ). More broadly, we have applied a similar 2 — a 

175 analysis to all ions in that report, at all temperatures and densities, and find the same 

176 agreement (20-26% confidence level). 

177 Lastly, it is of interest to investigate how our final suppression factor in Eq. [8] compares 
17B to our original, motivating, formulation of Eq. [8] for C IV, shown in Fig. [H There is 

179 generally goo d qualitative agreem ent. However, it is seen that the original formulation. 



based on the 



Badnell et al. 



( 1l993[ ) results, shows a somewhat stronger suppression effect 



181 up to X ^ 11. This is likely due to the more accurate treatment of the partial DR data of 



Badnell et al. 



( Il993[ ) entering the coUisional-radiative modeling, although some difference 



183 due to the coUisional-radiative modeling itself may also be present. This indicates that even 
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Fig. 2. — A comparison b etween th e 



present para meterized suppression factor and the 



colhsional radiative results of lSummersI (11974 fc 1979[ ) for a sample of ions and temperatures, 
as a function of density. 
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Fig. 3. — Agreement between the suppression curve of Eq. [9]and the 
data for all iron ions Fe''"'", g = 9 — 19. The upper panel shows the detailed level of agreement 
of the two end cases, Fe^"*" and Fe^^"*". The lower panel shows the 2 — a (95.4%) confidence 
as a function of charg e state; this means that 95.4% of all density data points in 



leve^ 

the ISummersI (11974 fc 19791 ) data, for the given charge and temperature, are within that 
percentage of the prediction from Eq. [HI The symbols denote different values of the scaled 
temperature 6 = T/q^. 
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Fig. 4. — A comparison between the final suppression factor of Eq. [8] (solid line), corre- 



sponding to a pseudo-Gaus sian profi 



and width w = 5.64548, the 



e with activat ion density Xa = 0.8314 (go = 40.284) 



SummersI ( 1l974 fc 19791 ) data points (solid diamonds), and the 



original formulation of Eq. |5] (dashed line), corresponding to a pseudo-Voigt profile with 
activation density Xa = 0.608, width w = 4.696, and mixture coefficient fi = 0.372. The 
temperature T = 1 x 10^ K is the same as in Fig. [H 
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184 collisional plasmas require coUisional-radiative modeling with state-of-the-art DR data. 

185 The stronger suppression tails-off at x > 11 as three-body recombination starts to become 

186 relevant and which, at even higher densities (not shown), causes the suppression factor to 

187 rise (since it is a ratio of effective recombination rate coefficients, i.e. includes three-body 

188 recombination.) 



2.3. Suppression Formula at Low Temperatures. 



The prec eding formulation was based on the suppression factor found by [Summers 



191 ( 1l974 fc 19791 ) for electron collisionally ionized plasmas, i.e., at higher temperatures, where 

192 DR is dominated by high-n resonances attached to a dipole-allowed core excited state. 

193 In photoionization equilibrium, however, the temperature at which a given ion forms is 

194 substantially smaller than that found in the electron collisional case. Due to the lower 

195 kinetic temperatures, DR occurs mainly through energetically low-lying autoionizing states, 

196 often via non-dipole core-excitations for which radiative stabilization is by the (outer) 

197 Rydberg electron. These states are not, in general, as susceptible to density suppression as 

198 their high-n counterparts, and so it may be necessary to modify the preceding suppression 

199 formulation. 

200 We first consider sequences with partially-occupied p-subshells in the ground state, 

201 which includes the B-like 2pi^Pi/2,3/2), C-like 2p^{^Po,i,2), 0-like 2p^(3Po,i,2), F-like 

202 2p5(2p3/2,i/2), Al-like 3p(2Pi/2,3/2), Si-likc 3p^CPo,i,2), S-likc 3p4(3Po,i,2), and Cl-like 

203 3p^(^P3/2,i/2) systems. For these sequences, there is fine-structure splitting in the ground 

204 state and a correspondingly small excitation energy, e^, giving dielectronic capture into high 

205 principal quantum numbers (because of the Rydberg relation g^/n^ < ctv)- Stabilization 

206 is via n ^ n' transitions and the recombined final state is built upon an excited parent. 

207 Ultimately, it is the strength of collisional coupling of this final state with the continuum 
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20B which determines whether recombination or ionization prevails. As the density increases, 

209 colhsional LTE extends further down the energy spectrum. It is difficult to give a general 

210 statement about the position of such final states relative to the ionization limit. So, 

211 we assume a worst case scenario, i.e., that such states are subject to suppression, and 

212 we use the preceding suppression formula. If density effects are found to be small in 

213 photoionized plasmas then this is likely sufficient. If they appear to be significant then 

214 a more detailed treatment based on coUisional-radiative modeling will be needed. Thus, 

215 for these systems, we retain the same suppression formula developed above, that is, 

216 S^{x, q, T) = S{x, Xaiq, T)) for = {5, 6, 8, 9, 13, 14, 16, 17}, and for all q and T. 

217 For the hydrogenic and the closed-shell He-like and Ne-like cases, on the other hand, 

218 the excitations proceed via an increase in core principal quantum number — Is —i- 2s or 

219 {2s, 2p} — {3s, 3p, 3d} — giving the dominant dielectronic capture into the low-n < 10 

220 resonances. Even following core radiative stabilization, these low-lying states are impervious 



221 to collisional reionization for the range of densities x < 10, anc 



thus we set S (x, q, T) = 1.0 



222 for = {1, 2, 10}. However, at densities x > 10, the ISummersI (11974 fc 19791 ) data for these 

223 three isoelectronic sequence show suppression factors that are fit well by the usual Eq. [HI so 

224 we do not modify S'^{x,q,T) for these cases. 

225 Lastly, we consider the intermediate isoelectronic sequences for which excitation 

226 arises from neither a fine-structure splitting of the ground state nor a change in principal 

227 quantum number of the core. These include the Li-like 2s — 2p, Be-like 2s^ — 2s2p, 

228 N-like 2s^2p^{^S) 2s2/(^P), Na-like 3s ^ 3p, Mg-like 3s^ 3s3p, and P-like 

229 3s^3p^{'^S) —7- 3s3p^(^P) cases up through the third row sequences. Any large low- 

230 temperature DR contribution arising from near threshold resonances is to low-lying states, 

231 for which suppression is negligible, i.e. the high-temperature suppression factor must be 

232 switched-off (S*^ — 1) at low-T. 
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233 To illustrate the general demarcation between low-T and high-T DR, we first consider 

234 DR of C IV, an overview of which is depicted in Fig. |3 The DR cross section, shown in the 

235 inset, is dominated by two features. The first is the n ^ oo accumulation of resonances 

236 at the 6 = 8 eV series limit — those whi ch can be treated in the usual high-T fashion 



237 ( iBurgess 



1965 



Burgess fc Summers 



19691 ) and are therefore susceptible to suppression 
23B according to our formulation above. However, there is a second strong contribution from 
239 the lowest accessible resonances just above the threshold electron energy, which, according 



240 to the Rydberg consideration S^/n^ ^ 63 = 0.6 Ryd, occur here f or n = 4. More genera. 



241 these low-lying states are typical of the low-lying DR spectrum (INussbaumer fc Storey 



1984 



ly, 



The ls^2p4/ resonances decay predominantly to the ls'^2s2p, ls^2p^ and ls^2s4/ 

243 states. These states lie well below the ionization limit and so are not susceptible to 

244 further reionization. Since there should be no density suppression then, we seek a modified 

245 suppression factor which tends toward unity (i.e., no suppression) at lower temperatures. 

246 In order to make a smooth transition from the high-T suppression factor S (x; q,T) 

247 given in Eq. [8], which is appropriate for the high-T peak region kT ^ kTmax = ^.e^/S, to the 

248 low-T region, where — )■ 1, we use the modified factor 



S'\x;q,T) 



[1- S {x; q, T)] exp 



lOkT ' 



(14) 



249 where eAr(g) = 8 eV for the particular case of C IV (A^ = 3 and q = 3). As seen in Fig. [5l the 

250 density- dependent effective DR rate coefficient, a'j^l^{ne,T), indeed s atisfies the requireme nt 



251 that the high-T peak is suppressed according to the formulation of 

252 whereas suppression is totally turned off for the lower-T peak. 



Badnell et al 



(119931) 



We have investigated the application of Eq. [14] for all ions that exhibit these same low-T 



^ We note that the C IV n = 4 res onance manifo^ 
near-threshold density-dependent effects (jPindzola et al. 



■ d has been the subject of further 
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Fig. 5. — DR of C IV. The inset shows the (zero-density) DR cross section convoluted with 
a 0.1 eV FWHM Gaussian. The spectrum is dominated by two features: the n = 4 DR 
resonance manifold below 1.0 eV and the n — )■ oo Rydberg resonances accumulating at the 
2s —J- 2p series limit 63(^0) ~ 8 eV. The main figure shows the effective DR rate coefficient 
for several densities. Our modified suppression formulation for x > 0, using Eqs. |8] and [T^ 
ensures that the high-T peak, corresponding to the n — )■ 00 Rydberg series of resonances, is 
suppressed but the low-T peak, corresponding to the n = 4 resonances, is not suppressed. 
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254 resonances features, namely, all isoelectronic sequences = {3, 4, 7, 11, 12, 15}, and we have 

255 found that the correct transitioning from suppression at the high-T-peak to no suppression 

256 at low-T is indeed satisfied, provided, of course, that the appropriate dipole-allowed 

257 excitation energy eAr(g) is employed. For efficient representation, the excitation energies 

258 along each isoelectronic sequence are parameterized by the expression 

5 



j=0 



q_y 
10/ 



(15) 



259 These parameters, which are determined by fitt ing the above expression to the available 



260 NIST excitation energies (IRalchenko et al. 



20n|), are listed in Table El 



261 We note that all isoelectronic sequences and ionization stages are now included in 

262 this prescription — our final comprehensive model for treating DR suppression, albeit 

263 in a simplified fashion. For those ions with fine-structure splitting in the ground state, 

264 we have eiy{q) ~ 0, so that S'^{x;q,T) = S{x;q,T). (We apply this generally also for 

265 Ar-like s equenc es and above (A^ > 18), based-on the density of states — see, for example. 



266 iBadnelll (j2006bl ).) For the closed-shell cases, on the other hand, we have eAr(g) — )■ oo. Thus, 

267 S'^{x;q,T) = 1 for hydrogenic and closed-shell systems, i.e., there is no suppression (for 
26B X < 10). Lastly, for the intermediate cases, the suppression factor is gradually increased 



269 toward unity at lower temperatures and begins to admit low-n DR resonances. 



270 3. Results 

271 The suppression factors derived here have been applied to the state-of-the-art total 

272 DR rate coefficients taken from the most recent DR database^ These modified data have 

273 been incorporated into version C13 of the plasma simulation code Cloudy, most recently 



http : / /amdpp . phys . strath . ac . uk/tamoc/DATA/ 
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Table 1. Modified A{N) coefficients from Eq. (fT2l) . 



Sequence A{Ny 



H-like 1 16 

He-like 2 18 

Li-like 3 66 

Be-like 4 66 

B-like 5 52 



^ These must be multiplied by 2.0 ii 9 < 2.5 x lO* K 



Table 2. Fitting coefficients for the excitation energies eAr(g) = Yl^j=oPN,j ^V. 
Numbers in square brackets denote powers of 10. 



Sequence 


N 


PN,0 


PN,1 


PN,2 


PJV,3 


PNA 


PN,5 


Li-like 


3 


1.963[-f0] 


2.030[-l-l] 


-9.710[-1] 


8.545[-l] 


1.355[-1] 


2.401 [-2] 


Bc-like 


4 


5.789[+0] 


3.408[+l] 


1.517[+0] 


-1.212[-f0] 


7.756[-l] 


-4.100[-3] 


N-like 


7 


1.137[-t-l] 


3.622[+l] 


7.084[-|-0] 


-5.168[-l-0] 


2.451[+0] 


-1.696[-1] 


Na-likc 


11 


2.248[-l-0] 


2.228[+l] 


-1.123[-l-0] 


9.027[-l] 


-3.860[-2] 


1.468[-2] 


Mg-likc 


12 


2.745[-l-0] 


1.919[-|-1] 


-5.432[-l] 


7.868[-l] 


-4.249[-2] 


1.357[-2] 


P-like 


15 


1.428[+0] 


3.908[+0] 


7.312[-1] 


-1.914[-f-0] 


1.051[-l-0] 


-8.992[-2] 


H-, He-, Ne-likc 


1,2,10 


oot 


0.0 


0.0 


0.0 


0.0 


0.0 


B-, C-, F-like 


5,6,8,9 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


A1-, Si-, S-, Cl-likc 


13,14,16,17 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 




> 18 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



t Reset to 0.0 for x > 10. 
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274 described by iFerland et al.l ( 120131 ) . Cloudy can do simulations of both photoionized and 

275 coUisionally ionized cases, and we show the effects of coUisional suppression on both. 

276 Figure [6] shows the ionization distribution of iron for the collisional ionization case. 

277 Figure [7| shows a similar calculation for photoionization equilibrium. Both show two 

278 hydrogen densities, 1 cm"'^, where collisional suppression of DR should be negligible, and 

279 1 0^" cm""^, where collisional suppression should greatly affect the rates for lower charges 

280 and temperatures. The upper panel shows the ionization fractions themselves, for these two 

281 densities, while the lower panel shows the ratio of the high to low density abundances. 

282 Cloudy's assumptions in com puting collisi o nal io nization equilibrium, as shown in 



283 Figure [6l have been described by 



Lykins et al. 



(120121 ). It is determined by the balance 



284 between collisional ionization from the ground state and recombination by radiative, 

285 dielectronic, and three body recombination to all levels of the recombined species. 

286 The photoionization case shown in Figure [7] depicts the Active Galact ic Nucleus 



287 spectral energy distribution (SED), described by 

288 of the ionization parameter 



Mathews fc Ferlandl (119871 ). function 



U 



riH c 



(16) 



289 where is the hydrogen-ionizing photon flux, Uh is the density of hydrogen, and c is the 

290 speed of light. There is only an indirect relationship between the gas kinetic temperature 

291 and the ionization of the gas in this case. Here, the level of ionization is determined by a 

292 balance between photoionization by the energetic continuum and the total recombination 

293 rate. 

294 The lower panels of Figs. [6] and [7] show that the amount that the ionization increases 

295 due to DR suppression can be large — the ratio can easily exceed 1 dex. Clearly, these 

296 results demonstrate that density effects on the ionization balance need to be considered 
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iron 




log Electron Temperature (K) 



Fig. 6. — Upper panel: coUisional ionization fractional abundance vs. electron temperature 
for all ionization stages of Fe. The solid curves correspond to a density of 1 cm~^ and the 
dashed curves correspond to a density of 10^° cm~'^. From left to right, the curves range 
from Fe I to Fe XXVII. Lower panel: ratio of the calculated fractional abundances for the 
two densities. 
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iron 



l.OOF 




-4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



logU 

Fig. 7. — Upper panel: photoionization fractional abundance vs. the ionization parameter 
U for all ionization stages of Fe. The solid curves correspond to a density of 1 cm~'^ and 
the dashed curves correspond to a density of 10^*^ cm"^. From left to right, the curves range 
from Fe I to Fe XXVII. Lower panel: ratio of the calculated fractional abundances for the 
two densities. 
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297 more precisely. 



298 4. Conclusion 

299 We have investigated the effects of finite densities on the effective DR rate coefficients 

300 by developing a sup pression factor model, which was motivated by the early work of 



Badnell et al. 



( I1993I ) for C IV and extended to all other ions using physicall y-motivate d 



302 scaling consi derations, and more precise fitting of collisional-radiative data fjSummers 



303 I1974 fc 19791 ) ■ Accurate zero-density DR rate coefficients were then multiplied by this 

304 suppression factor and introduced into Cloudy to study the finite-density effects on 

305 computed ionization balances of both coUisionally ionized and photoionized plasmas. It 

306 is found that the difference in ionization balance between the near-zero and finite-density 

307 cases is substantial, and thus there is sufficient justification for further studies of collisional 

308 suppression from generalized collisional-radiative calculations. This is expected to impact 

309 the predictions of the ionization balance in denser cosmic gases such as those found in nova 

310 and supernova shells, accretion disks, and the broad emission line regions in active galactic 

311 nuclei. 

312 The present results are intended to be preliminary, and to demonstrate the importance 

313 of density effects on dielectronic recombination in astrophysical plasmas. Given the 

314 approximations adopted, we suggest that their incorporation into models (e.g., via Cloudy) 

315 be used with a little caution. For example, one might run models with and without the 

316 effects of suppression at finite density, especially in modeling higher density plasmas (e.g.. 



317 the broad emission 



i ne reg ion in quasars). Nevertheless, it is nearly half a century since 



318 iBurgess k. Summers! (119691 ) demonstrated significant density effects on DR, and it is time 

319 that some representation exists within astrophysical modeling codes to assess its impact on 

320 the much more rigorous demands made by modern day modeling, especially given its routine 
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321 incorporation by magnetic fusion plasma modeling codes. In the longer term, we intend 

322 to present results based on detailed coUisional-radiative calculations using state-of-the-art 

323 state-specific DR rate coefficients. 



324 5. Acknowledgments 

325 DN, TWG, and KTK acknowledge support by NASA (NNX11AF32G). GJF 

326 acknowledges support by NSF (1108928; and 1109061), NASA (lO-ATPlO-0053, 10- 

327 ADAPlO-0073, and NNX12AH73G), and STScl (HST-AR-12125.01, GO-12560, and 

32B HST-GO-12309). UK undergraduates Mitchell Martin and Terry Yun assisted in coding the 

329 DR routines used here. NRB acknowledges support by STFC (ST/J000892/1). 



- 26 - 

330 REFERENCES 

331 Badnell, N. R. 2006a, ApJS, 167, 334 

332 Badnell, N. R. 2006b, ApJ, 651, L73 

333 Badnell, N. R., O'MuUane, M. G., Summers, H. R, et al. 2003, A&A, 406, 1151 

334 Badnell, N. R., Pindzola, M. S., Dickson, W. J., et al. 1993, ApJ, 407, L91 

335 Bates, D. R., Kingston, A. E., & McWhirter, R. W. R 1962, Proc. R. Soc. Lond. A, 267, 

336 297 

337 Burgess, A. 1965, ApJ, 141, 1588 

338 Burgess, A. & Summers, H. R 1969, ApJ, 157, 1007 

339 Davidson, K. 1975, ApJ, 195, 285 

340 Ferland, G. J., Porter, R. L., van Hoof, P. A. M., et al. 2013, ArXiv e-prints 

341 Jordan, C. 1969, MNRAS, 142, 501 

342 Lykins, M., Ferland, G. J., Porter, R. L., et al. 2012, MNRAS, submitted 

343 Mathews, W. G. & Ferland, G. J. 1987, ApJ, 323, 456 

344 Nussbaumer, H. & Storey, P. J. 1984, A&AS, 56, 293 

345 Osterbrock, D. E. & Ferland, G. J. 2006, Astrophysics of Gaseous Nebulae and Active 

346 Galactic Nuclei, 2nd edn. (Sausalito, CA: University Science Books) 

347 Pindzola, M. S., Loch, S. D., & Robicheaux, F. 2011, Phys. Rev. A, 83, 042705 

348 Ralchenko, Y., Kramida, A. E., Reader, J., & NIST ASD Team. 2011, National Institute of 



Standards and Technology, http://physics.nist.gov/asd 



350 Summers, H. P. 1972, MNRAS, 158, 255 

351 Summers, H. P. 1974, MNRAS, 169, 633 

352 Summers, H. P. 1974 & 1979, Appleton Laboratory Internal Memorandum IM367 & 

353 re-issued with improvements as AL-R-5 

354 Summers, H. P. & Hooper, M. B. 1983, Plasma Physics, 25, 1311 

355 Verner, D. A. & Yakovlev, D. G. 1995, A&AS, 109, 125 

356 Voronov, G. S. 1997, Atomic Data and Nuclear Data Tables, 65, 1 

357 Wertheim, G. K., Butler, M. A., West, K. W., & Buchanan, D. N. E. 1974, Rev. Sci. 

358 lustrum., 45, 1369 



This manuscript was prepared with the A AS lATgX macros v5.2. 



